THERMO-ELASTICITY FOR ANISOTROPIC MEDIA 
IN HIGHER DIMENSIONS 



JENS WIRTH 



Abstract. In this note we develop tools to study the Cauchy problem for the sys- 
tem of thermo-elasticity in higher dimensions. The theory is developed for general 
homogeneous anisotropic media under non-degeneracy conditions. 

For degenerate cases a method of treatment is sketched and for the cases of cubic 
media and hexagonal media detailed studies are provided. 



1. Introduction 

While isotropic thermo-elasticity is a well-known and well-established subject (see, 
e.g., the book of Jiang Song-Racke [TS] and references therein) only very few results are 
available for the case of anisotropic media. Among them are the theses of Borkenstein 
[1] for cubic media and Doll [3] for the case of rhombic media together with the authors 
treatments [16], [22], all in two space dimensions. 

In this paper the system of anisotropic thermo-elasticity in three (and more) dimen- 
sions, i.e., 

U tt + A(D)U + jV9 = 0, (1.1a) 
0t-«A0 + 7V-E/'t = O (1.1b) 

for the elastic displacement U(t, •) : W 1 —> M. n and temperature difference 6(t, •) : M. n — > R 



to the equilibrium state, will be considered. The system (1.1) couples the hyperbolic 
elasticity equation with the parabolic heat equation. The operator ^4(D) describes the 
elastic properties of the underlying medium, while k denotes its thermal conductibility. 
The constant 7 describes the thermo-elastic coupling. Basic assumptions of our theory 
are k > 0, r y 2 > together with 

• A(£) = \^\ 2 A(r]), 7] = £/|£|, is a 2-homogeneous matrix-valued symbol; 

• A : S n_1 -> R nxn is a real-analytic function of 77 G S rt_1 , n > 3; 

• A(rj) = A*(tj) > is self-adjoint and positive. 

In general we can not assume that A{rj) is non-degenerate in the sense that # spec A{rj) = 
n for all 77 £ S™ _1 (as done for the two-dimensional case in [16 1). All basic examples 
show degeneracies in dimensions n > 3. 

Example 1.1. Isotropic media 

A{rj)= \il +(A + n)r]®r} (1.2) 



Key words and phrases, thermo-elasticity, a-priori estimates, anisotropic media, degenerate hyper- 
bolic problems. 
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with Lame constants A and \x. The matrix A(rj) is positive as long as /j > and A > —2fi. 
The eigenvectors of A(rj) are multiples of r\ and n 1 - and thus invariant under rotations 
of frequency space. 



Example 1.2. Cubic media 

({t - n)rH + n 

Mm- ( A + ^2 



(A + ^mm 



(A + /U)77i77n \ 



(1.3) 



V (^ + A t ) r ?l ? 7n ••• ••• (r- ii)rft + fij 

described by parameters A, \x and r. Later we will describe the assumptions made on 
these parameters and the resulting spectral properties of the matrix function A(rj) more 
precise. In the case of three space dimensions, the matrix A(rj) is positive if and only if 
fx > 0, r > together with — 2/j — t/2 < A < r. In three space dimensions this will be 
one of our main examples. 



Example 1.3. We can replace the constant r on the diagonal by t\, . . . , r n in ( 1.3). This 



yields so-called rhombic media. The behaviour of rhombic media is close to that of 
cubic media if the parameters are of similar size, in general there will appear exceptional 
situations. See, e.g., |22j or [21] for a discussion of this effect in two space dimensions. 



Example 1.4. Hexagonal media are another particularly interesting case for three space 
dimensions. Since we want to come back to them later on we introduce the corresponding 
operator. It is given by 

A(rj) = B(7]) T CB(r]), (1.4) 

where C contains the 5 structure constants t\, t%, Ai, A2 and \i and D(r/) is of a particular 
form, 
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(1.5) 



/ 



Even the first (non-trivial anisotropic) example, the case of cubic media in three space 
dimensions, has degenerate directions in which A(rj) has double eigenvalues. Later on 
we will analyse this example in detail. 

Definition 1. We call a direction r\ S S ra_1 (elastically) non-degenerate if 

# spec ^4(77) = n (1-6) 

holds true for this direction 77. 

The set of non-degenerate directions is an open subset of S™" 1 . For non-degenerate 
directions the treatment of |16j transfers almost immediately and gives a representation 
of solutions. We will sketch the results in Section [2] In Section [3] we consider special 
degenerate directions and discuss the examples of cubic and hexagonal media. Dispersive 
estimates for solutions are given in Section |4| In the neighbourhood of degenerate 
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directions they are essentially based on estimates developed by Liess [8], [10] for the 
treatment of anisotropic acoustic equations. 

2. Treatment of non-degenerate directions 

For the following we consider a simply connected open subset U of S n_1 , where the 
symbol A(rj) has n distinct (and real) eigenvalues. We denote these eigenvalues in 
ascending order as 

< x x (?7) < x 2 (r?) < ••• < x n (rj). (2.1) 

By analytic perturbation theory, see [6], we know that these eigenvalues are real-analytic 
and that we find corresponding normalised eigenvectors 

r 1 (r ? ),...,r n (7 7 ) eC^iU^ 1 - 1 ) 
depending analytically on rj £ U. Collecting them in the unitary matrix 

M (v) = { r i(r])\ r 2(v)\ ■ ■ ■ \ r n(v)), 
M*(rj)M(rj) = I = M(r])M*(rj), 

we can diagonalise the matrix A{rj) 

A(rj)M(r)) = M(r/)P(r/), 

V(rf) = diag (xi(rj), x 2 (v),---, ^nO?))- 

In our treatment we will not make use of analyticity directly, instead our use of pertur- 
bation theory will be based on [5] und [23] and uses only smooth dependence. This will 
be of interest for generalisations later on. Therefore, whenever we use analyticity, we 
will explicitly state that. 

We use M(rj) to reduce the thermo-elastic system to a system of first order. For this 
we denote by U and 9 the partial Fourier transforms of U and 9 with respect to the 
spatial variables and consider 

'(Dt + 2?V3(£))M*fa)#\ 
(T> t -V^))M*{T,)u\ G 



(2.2) 

(2.3) 
(2.4) 

(2.5) 
(2.6) 



V 



i2n+l 



(2.7) 



as usual D( = —\dt and r\ = Then V satisfies a first order system of ordinary differ- 

ential equations, which has an apparantly simple structure. Straightforward calculation 
shows that 



D t V = B(£)V 



(2. 



holds true with coefficient matrix 

"2(0 



ljai(e) JoaCO 



-"2(0 



WO 



vya 2 (0 

1701(6 
1702(6 



(2.9) 
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where Wj-(f) = y^I?) G C°°(W,IR+) and 

a 3 {t) = r M<- ( 2 -10) 
Following the conventions of [16J we denote these functions aj{£) as the coupling func- 
tions of the thermo-elastic system associated to the elastic operator A(D). They play a 
prominent role for the description of the time- asymptotic behaviour of solutions. This 
reflects the fact that they couple the homogeneous first order entries in B(£) with the 
second order lower right corner entry. Note, that 

n 

Zeros of the coupling functions are of particular importance. Following |16|. Def. 1] we 
define: 

Definition 2. A non-degenerate direction r\ G § n ~ 1 is called 

• hyperbolic if one of the coupling functions vanishes; more precisely, it is called 
hyperbolic with respect to the eigenvalue Xj{rf) if 07(77) = 0; 

• parabolic if all coupling functions are non-zero. 

In the anisotropic case the set of hyperbolic directions is (genericallj^]) a lower dimen- 
sional subset of S n . In order to decide whether a direction is hyperbolic or parabolic 
we can employ the following proposition. We denote for a matrix A and a vector r\ by 

Z(A, rf) = span{ A k rj\ k = 0,1, .. .} (2.12) 

the corresponding cyclic subspace, i.e. the span of the trajectory of r\ under the action 
of the matrix A. 

Proposition 2.1. The following statements are equivalent: 

(1) The cyclic subspace of n has dimension n — k, i.e., &vca.Z(A(rj), rj) = n — k. 

(2) Exactly k of the coupling functions vanish in rj. 

Hence, a non- degenerate direction r\ G is parabolic if and only if Z{A(n),n) = W l 

and therefore 

det (n\A(n)n\ ■ ■ ■ \A n - l {n)r,) / 0. (2.13) 
Proof. If we represent 77 in the eigenbasis of A{rf) we obtain 

n = ai(rj)ri(rj) H h a n (n)r n (n) (2.14) 

and therefore 

A l (n)n = ^(77)01(77^1(77) + • • • + x%)a n {n)r n (ri). (2.15) 

If k of the coupling functions vanish, then A n ~ k {7])n must be in the span of the ^(77)77 
with t = 0,1,. ..,71 — k — 1 and thus the cyclic subspace is at most of dimension 
n — k. On the other hand, the first n — k vectors in the trajectory are linearly in- 
dependent since the corresponding matrix in the basis representation with respect to 



^"If not, by analyticity it follows that one coupling function vanishes on U and the system is therefore 
decoupled. This case is reduced to the study of the lower dimensional blocks, one is a hyperbolic system 
the other one a thermo-elastic system of lower dimension. This is, e.g., the case for hexagonal media, 



see Section 3.4 
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ai (77)7*1 (77), . . . , a n (rj)r n (rj) is just the van der Monde matrix associated to the eigenval- 
ues of A(r]) for non-vanishing coupling functions and therefore regular. □ 

2.1. On the characteristic polynomial of the full symbol. At first we collect 
some of the spectral properties of the matrix -B(£) which are directly related to the 
characteristic polynomial of B(£). 

Proposition 2.2. The following identies hold true: 

tiB(0=i^\ 2 , (2.16) 
det 5(f) = ik|£| 2 det (2.17) 

n n 

det(i/ - B(0) = (u- i^| 2 ) H(v 2 - " W £ a "M) II^ 2 - **(*))■ (2-18) 

3=1 j=l k^j 

Furthermore, the matrix B{£) has a purely real eigenvalue for £ 7^ if and only if the 
direction 77 = £/|£| is hyperbolic. If it is j -hyperbolic, then ±Wj(f) G spec-B(f). 

The proof of the last fact is fairly straightforward and consists of separating real and 
imaginary parts of the characteristic polynomial. Note that for all parabolic directions 
we can divide the characteristic polynomial by i / Y[j( 1 ' 2 ~ to obtain 



i K |C| 2 o -A «K0 , , 

1 = — ^- + 7 2 V 9 J ' ■ 2.19 

This formulation allows to consider the neighbourhoods of hyperbolic directions. Assume 
for this that the set of hyperbolic directions with respect to Xj (77) 

Mj = {77 G U I aj{rj) = rj(ri) ■ 77 = 0} (2.20) 

is a regular submanifold of U. If we consider the corresponding hyperbolic eigenvalues 
°f B(£) in a neighbourhood of Mj, i.e. the eigenvalues which satisfy 

lim uf(0 = ±^(0 (2.21) 



for fixed |f |, equation (2.19) gives a precise description of the behaviour of the imagi- 
nary part of these eigenvalues. The proof is a straightforward generalisation from |16l 
Prop. 2.2]. 

Proposition 2.3. The non-tangential limit 

* J' 7* = 1 ^ " ^ E J <() m = T (2.22) 

^+ A/, z/± (f ) 2 _ Xj . (£) Wj . (£) ^ Xj (£) - Xfc (£) 

exists and is non-zero for all f 7^ 0. Furthermore, 

lim 2 ^ = - . . .^y 1 > 0. (2.23) 

u-^i a 2 (7?) 2w,(7 / )(C? + |£| 2 D?) 
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Figure 1. Non-zero eigenvalues of -Bi(^) for parabolic directions. 



2.2. Asymptotic expansion of the eigenvalues as |£| — > 0. We decompose B(£) 
into homogeneous components B(£) = i?i (£ ) + i?2 (£) of degree 1 and 2, respectively. For 
sufficiently small |£| we expect the eigenvalues of B(£) to be close to the eigenvalues of 
B\(£). For parabolic directions the (non-zero) eigenvalues of Bx(rj) can be determined 
from the equation 

1 A ajfo) 



r 



Z> 2 



(2.24) 



which follows directly from (|2.19|) with k 
Figure 



n 



0. It can be solved (e.g. graphically, see 
= 3) to obtain the distinct eigenvalues 0, ±£1(77), . . . , ±z> n (r?) ordered as 

< wifa) < ^i(ry) < uj 2 (v) < Mv) <•••< u n {v) < v n {rj). (2.25) 

For hyperbolic directions a similar result holds true. In the case of hyperbolic direc- 
tions w.r.to Xj{rj) eigenvalues move to ujj(rj). According to the choice of the coupling 
constant 7 different cases occur: 

(1) if -y is large then Uj(rj) = (jJj{rj), the other inequalities are unchanged; 

(2) if -j is small then £^-_]_ (77) = ojj(rj) and the other inequalities remain true. 

The critical threshold between these two cases is 

,2, 



1 



t 2 ^ xj(v) - *k(n) 



where -61(77) has the double eigenvalue l>j-\{rj) 
tions from 1161 we define: 



(2.26) 

Vj(jl)- Following the conven- 



Definition 3. We denote a hyperbolic direction w.r.to Xj{rf) as 7-degenerate if (2.26) 
holds true. 

For the following treatment we exclude 7-degenerate hyperbolic directions and assume 



instead that for all hyperbolic directions in U. condition (2.26) is not satisfied for the 



corresponding index j. Then the following statement is apparent. 
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Proposition 2.4. Let rj be not 7- degenerate. Then the matrix B\(rj) has 2n + 1 distinct 
real eigenvalues 0, ±17, . . . , ±z/ n for all r\ £U. 



Proposition |2.4| allows to apply the standard diagonalisation scheme (see Sec. 2.1]) 
to = B\ (£) + £>2 (£) as £ — ^ 0. Hence, eigenvalues, eigenprojections and all their 

derivatives have full asymptotic expansions as £ — >• 0. The proof is almost identical to 
that from |16} Prop. 2.5] and is omitted. 

Proposition 2.5. For all not ^-degenerate directions 77 = £/|£| G £/te eigenvalues and 
eigenprojections of B{£) have full asymptotic expansions as £ — > 0. T/ie mam terms are 
given by 

MO = ^\ 2 MV) + Om 3 ) (2.27a) 
if (0 = ±|^fa) + H£| 2 ^) + Odei 3 ) (2.27b) 



with 



and 



-l 



-1 



Furthermore, bj{rf) = if and only if n is hyperbolic with respect to the eigenvalue 
Remark 2.1. Note, that tr_B(£) = ift|£| 2 implies 

n 

b (r ] )+2j2b j (v) = l- (2.29) 



Recall that by Proposition 2.2 eigenvalues of B(£) can only be real along hyperbolic 
directions (and then they are exactly the 'trivial' real eigenvalues). In combination with 
the fact that eigenvalues of B(£) are continuous in £ we obtain: 

Corollary 2.6. For all parabolic directions r] = £/|£| 6 hi we have Imi/- (£) > 0. TTie 
same is £rue as /ong as 77 is not hyperbolic w.r.to Xj(rj). 

2.3. Asymptotic expansion of the eigenvalues as |£| — > 00. In this case the two- 
step procedure developed in |5, Sec. 2.2], [TBI Prop. 2.6] applies in analogy. Essential 
assumption is the non-degeneracy of A(jj). We omit the proof and cite the corresponding 
result only. 

Proposition 2.7. For all non- degenerate directions the eigenvalues and eigenprojections 
of the matrix B(£) have full asymptotic expansions as |£| — > 00. The first terms are given 
by 

^o(O = ^l£| 2 -- + 0(ir 1 ), (2.30a) 

K 

uf(0 = ±\t\ Uj { V ) + ^-a)(r,) + Od^- 1 ). (2.30b) 
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Remark 2.2. Despite the fact that we used the same notation for the eigenvalues as 
£ — > and |£| — > oo, we do not claim that they are indeed the same functions of £. 
This is only true for hyperbolic eigenvalues near hyperbolic directions, in general there 
might be multiplicities in between and there might be no consistent notation for these 
functions. 

Corollary 2.8. For all parabolic directions n = £/|£| the eigenvalues of B(^) satisfy 
lmv(r]) > C v > for |£| > c. The same is true for parabolic eigenvalues in hyperbolic 
directions. 

Remark 2.3. In particular, we see by the asymptotic expansions that the eigenvalues 
of are simple for large and also for small values of |£|. Furthermore, we see that 
the hyperbolic eigenvalues are always separeted (i.e. if multiplicities occur in hyperbolic 
directions, they involve only parabolic eigenvalues). 



2.4. Behaviour of the imaginary part. The asymptotic expansions of Propositions 2.5 



and 2.7 allow to draw conclusions for the behaviour of the imaginary part. We collect 
them for later use. The first result is apparent. 

Proposition 2.9. On any compact set of parabolic directions we have the uniform esti- 
mates 

Imuf ] (C)>C e forall\£\>e, (2.31) 

Wf \0 ~ b 3 ( V )\d\ 2 for all |C| < e (2.32) 

for all eigenvalues of B{^) and arbitrary e > 0. 

The next statement is concerned with a tubular neighbourhood of a compact subset 
of a regular submanifold Ma of hyperbolic eigenvalues w.r.to >Cj(rj). It is only of inter- 
est how the corresponding hyperbolic eigenvalues ^f(0 behave, the others still satisfy 
Proposition |2.9[ 



Proposition 2.10. Uniformly on any tubular neighbourhood of a compact subset of Mj 
of non-j- 
estimates 



of non-j- degenerate directions the corresponding hyperbolic eigenvalues ir- (£) satisfy the 



Imi/±(0 ~ a)^) for all |£| > e, (2.33) 

lm^(0 ~ bj(r,m 2 for all |£| < e. (2.34) 

Proof By Proposition |2.3| we know that 

hnif(£) = aj($)K{$ (2.35) 

for some function K(£). Our aim is to estimate K(£). The left hand of this formula has 
a full asymptotic expansion as |£| — > and |£| — > oo. Therefore, also the right hand side 
has one and it follows that 

m) = t~ + °(l£r 1 )> ICI^oo, (2.36a) 



m) = ^\ 2 ^M + o(ier 3 ), ic i o. (2.36b) 

a z Arj) 
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Thus, the desired estimate follows by a compactness argument as soon as we have a uni 
form lower/upper bound for bj (r/)/a| (17) • The representation of bj(rj) in Proposition 
in combination with ( |2.24 ) implies 



2.5 



hm — — — = hm 7 (i/- +Xj(?7)) 7 xn r~^9 



=2 7 2 x J (f?)fl- 7 2 V °* ( ^ , (2.37) 

which is clearly bounded and (uniformly) positive on any compact subset of Mj (where 
we have to use that 77 G Mj is not 7-degenerate) . □ 

2.5. Conclusions. We will draw several conclusions what we have obtained so far and 
what we still have to consider in the remaining part of this treatise. 

2.5.1. Cubic media in 3D. If we consider the special case of cubic media in three space 
dimensions degenerate directions are given by 77 = (771, 772, 773 ) T with fjf = fj 2 = ffe (eight 
directions, corresponding to the corners of a cube) or fjf = 1 for some i (six directions, 
corresponding to its faces). This can be calculated directly, corresponding eigenspaces 
are span{77} and fj 1 - = {£ G W 1 \ 77 • £ = 0}, or concluded by the cubic symmetr^ of A(£) 
in this particular case. See Figure [2j 



To obtain the hyperbolic directions we apply Proposition 2.1 and look for the action 
of r] under ^(77). We obtain that 

(1) a direction 77 is hyperbolic if and only if 

det(7 ? |^(77)77|A 2 (7?)7?) = (r - A - 2 / x) 3 7 ?1 7/ 2 77 3 (77 2 - t ?2 2 )(t / 2 - t^t/ 2 - t/ 2 ) = 0, (2.38) 

thus the set of hyperbolic directions is the union of nine great circles on S 2 ; 

(2) 7711^4(77)77 for all 26 intersection points of these great circles, 14 of them are ex- 
cluded as being degenerate. 

Except for these 14 points on S 2 we obtained an almost complete description of the 
spectrum of B(£). We know full asymptotic expansions of eigenvalues for small and 
large frequencies |£|, estimates for the imaginary part of them and similar statements 
for eigenprojections. This information allows to draw conclusions on the large time 
behaviour of solutions, e.g. energy and dispersive estimates. This can be done similar 
to the treatment of [16] . see Section |4j The remaining degenerate directions appear in 
two types, which can be interchanged by the action of the symmetry group. The study 
of these degenerate directions is what is left open so far and will be the main point of 
Section |3j 



2 A(£) is invariant under the hexaeder group, i.e. the symmetry group of a cube. Thus, eigenspaces 
must be transferred in an appropriate way, which implies that symmetries of order 3 or 4 can only be 
realised by higher dimensional eigenspaces. 
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Figure 2. Degenerate points for cubic media correspond to symmetries 
of a cube. Corner points A are conic singularities, midpoints of faces 
B uniplanar singularities of spec A(rj). The midpoints of edges C are 
non-degenerate, but hyperbolic with respect to two different eigenvalues. 



2.5.2. Isotropic media. If we consider the special case of isotropic media, A{rj) = fil + 
(A + fi)r] <g) rj, we see that spec A(rj) = {fi, A + 2/j,} and corresponding eigenspaces are 
span {77} (corresponding to A + fi) and rj 1 - (corresponding to n). All directions are 
(elastically) degenerate. However, we still find locally smooth systems of eigenvectors. 
All directions are hyperbolic and the hyperbolic eigenvalue \x has multiplicity n — 1. 
Therefore the system D t V = B(£)V decouples into a diagonal part of size 2n — 2 and a 
full 3x3 block and is given after a rearrangement of the entries as 



B(t,0 



\ 



vX+2^|£| 



V 



-vX+2^1 
-m 



(2.39) 



/ 



This block structure corresponds to the Helmholtz decomposition of vector fields applied 
to the elastic displacement. If V ■ U(t, ■) = the lower block cancels and we obtain wave 
equations with speed JJi for the components of U. Otherwise, if we cancel the upper 
block we obtain the 3x3 system corresponding to one-dimensional thermo-elasticity 
with its well-known properties. 
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2.5.3. One-dimensional thermo- elasticity. For completeness we mention some results on 
the one-dimensional system 

u tt - t 2 u xx + j9 x = 0, (2.40a) 

6t - k6 xx + ju tx = 0. (2.40b) 

We assume 7, k, t > 0. Following our strategy we can rewrite this problem as first order 
system. The corresponding symbol B(£) is given by 

/ r£ ij£ \ 

B(0 = -t£ i 7 £ • (2.41) 

Its eigen values satisfy asymptotic expansions for £ — > and £ — > ±00. Propositions 2.5 
apply with P 1 * 1 = ±\/t 2 + 7 2 and 



and 



2.7 



^ 2 b _1 7 2 

r 2 + 7 2 ' 2 r 2 + 7 2 ' 



(2.42) 



Therefore, by Proposition 2.5 



M0 = i^^e + O{e), (2.43a) 
t z + 7 Z 

= iy^+T^ + i „, ^ 2 , g 2 + C(e 3 ), (2.43b) 
2(r z + Y) 



as £ — >• and by Proposition |2.7| 

MO = i< 2 -i- + 0(r 1 ), (2.43c) 

ft 

z, 1 ± (0 = ±re + i^ + O(r 1 ), (2.43d) 

as £ 00. The essential information for large time estimates is given by the behaviour 
of the imaginary part. It follows that Imz/(£) > C t for |£| > e for certain constants and 

W0~-£p^ a . W(0~2pq^' (2 - 44) 

2.5.4. Hexagonal media in 3D. For hexagonal media in three space dimensions the sit- 
uation is (surprisingly) simpler than for cubic media. The elastic operator defined by 
(1.4)-(1.5) is invariant under rotation around the £3-axis (taking into account a corre- 
sponding rotation of the reference frame for vectors) and therefore it suffices to under- 
stand its cross sections in the x\—xi plane. We will sketch some of the properties of the 
corresponding symbol -A (77). 

Following Proposition |2. 1| we obtain 

(1) that 

det(r/| A{n)n\ A 2 (77)77) = 0, (2.45) 
such that all directions 77 £ S 2 are hyperbolic. The corresponding eigenspace is 
(generically) given by multiples of (772, —771, 0) such that the hyperbolic eigenvalue 
is 

"^(vf + nD + wl (2-46) 
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(2) It remains to look for directions with two hyperbolic eigenvalues. They satisfy 
7/||^4(7/). This is true, if 773 = or if 7/1 = 772 = or if 

A 2 + 2/j 



4 



n 



(2 47) 

2A2 + 4/J + T1-T2' v ' ' 

provided the latter expression is non-negative. Except in the limiting case t\ = 
A2 + 2/j, the coupling functions vanish to first order along the corresponding 
circle. If t\ = r 2 = A2 + 2/j all directions are hyperbolic with two hyperbolic 
eigenvalues and if t\ = A2 + 2/j 7^ T2 coupling functions vanish to third order. 
(3) The matrices A(rj) are invariant under rotation. Introducing spherical coordi- 
nates on S 2 



cos 4> cos ip + 



sin <f> cos ifi + 



sin tp 



(2.48) 



and using a corresponding (moving) basis for vectors given by 



±1 



V 7 ?! + ^2 




( 



X 



'/I 

'/2 

>7i' 



(2.49) 



(sign chosen to make them smoothly dependent on 7/ / ±(0,0, 1) T ) decom- 
poses A(rj) into (l,2)-block-diagonal structure (independent of the co-ordinate 
<p). The scalar block corresponds to the eigenvalue (2.46), while the 2x2 
block has trace \x-\-t\ cos 2 ip + T2 sin 2 ip and determinant /xn cos 4 ip 



T\T2— 2A2 



jUT2 sin 4 -0 + 



sin 



2^. 



If (n — /j)(t2 — /x) 7^ 0, the 2x2 block has distinct eigenvalues for all -0 and therefore 



the only degenerate directions are directions where this block has (2.46) as one of its 
eigenvalues. This happens if and only if the right hand side of (2.47) is non-negative 
and on the circle defined by that equation. 

Thus, the previously developed theory is applicable for all directions except the de- 
generate ones 7/1 = 7/2 = or (2.47). The always existent hyperbolic eigenvalue (2.46) 



leads to a decoupling of the thermo-elastic system into two scalar blocks and a (at least 
formally) 2D thermo-elastic system. 

Due to rotational invariance, it suffices to treat the cut 771 = for handling of degen- 
erate directions. This will be sketched later. 



3. Some special degenerate directions 

We want to study neighbourhoods of degenerate directions for some particular cases. 
To study degenerate directions in full generality is beyond the scope of this paper. We 
relate our approach to the type of singularity of the corresponding Fresnel surface 

S = {£ G R n I 1 G spec A(£)}- (3.1) 

This surface is in general n-sheeted and for all non-degenerate directions these sheets 
are given by 

Sj = {£ G M n non-deg. | = 1} = {wj 1 (77)77 1 77 G S n_1 non-degenerate}, (3.2) 
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Figure 3. A cut through the Fresnel surfaces for examples of a cubic 
and a hexagonal medium. The material parameter are A = 1, r = 4 and 
\i = 1 for the picture on the left (cubic) and Ai = 1, A2 = g, t\ = 4, 
T2 = 1 and fj, = 3 for the picture on the right (hexagonal). 



while in degenerate points the surface is self-intersecting. For the importance of these 
surfaces in elasticity theory and some interesting properties of them we refer to Duff [4J 
or the investigations from Musgrave p2], |13| and Miller- Musgrave [14j . 

We remark only one of the general properties of S here. If A(£) is polynomial in £ 
then the surface S is algebraic of degree 2n and therefore any straight line intersecting 
S has at most 2n intersection points with S. In particular, if the inner sheet S n does 
not touch any of the the outer sheets, it has to be strictly convex. 



3.1. General strategy. If we investigate isolated degenerate directions or regular man- 
ifolds of degenerate directions of codimension greater than one we are faced with two 
major obstacles. Generically, eigenvectors of A (77) can not be chosen continuously in 
a neighbourhood of the degenerate direction and therefore a reformulation as system 
of first order as in (2.7) is problematic. This problem is related to higher-dimensional 
perturbation theory of matrices. It is well-known that in the one-dimensional situation 
eigenspaces are continuous (see, e.g., the book of Kato, [6]) and it can be resolved by 
introducing polar co-ordinates / normal co-ordinates around the degenerate directions 
and a system related to (2.7) can be formulated on a corresponding blown-up space (see, 
e.g., (3.10) below). Second obstacle are the multiplicities itself. Eigenvalues and eigen- 
vectors of the constructed system of first order do not possess asymptotic expansions in 
powers of |£| as |£| tends to or 00. However, especially in the three-dimensional setting 
we can write full asymptotic expansions in the distance to the degeneracy uniform in all 
remaining co-ordinates. 
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o 0.25 0.5 0.75 l -1 -0.75 -0.5 -0.25 -0 0.25 0.5 0.75 1 



Figure 4. Cuts of the Fresnel surface for cubic media; on the left hand 
side in the plane 773 = 0, on the right for 7/1 = 772. The parameters are 
chosen as r = 8, A = 2 and fi = 2. 



We will discuss the application of this general strategy in detail for conic singularities 
of the Fresnel surface appearing for the case of cubic media and give the corresponding 
results for uniplanar singularities afterwards. Finally we will consider hexagonal media 
and show that they are much simpler in their analytical structure. 



3.2. Cubic media, conic singularities. The Fresnel surface for cubic media has eight 
conic singularities which are related by the symmetries of the underlying medium. It 
suffices to consider one of them and we choose 77 = ^(1, 1, l) r £ §> 2 - Near this direction 

we introduce polar co-ordinates (e, 4>) on the sphere § 2 by 



cos <fi + e—= —1 sin (j). (3.3) 





They allow to blow up the singularity by looking at [0, 00) x S 1 instead of M 2 as local 
model of § 2 . In order to simplify notation, we apply a diagonaliser M of A(fj) to our 
matrices. For this we choose the unitary matrix 




M = —= V2 -1 -V3 (3.4) 



(corresponding to the vectors chosen already in ( |3.3| )). The matrix M 1 A(ri)M has a 
full asymptotic expansion as e — > and can be written as 



M- l A(e,(j))M = Ao + eAii^ + Oie 2 ), e^O 



(3.5) 
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with matrices 



,r + 2A + 4/i r + fj-A r + ^-A v , , 
= diag , , , (3.6a) 



. . .. 2r - /j + A 
^i(^) = £ 



3 3 

cos fli> sm ( 






+ — — '- — — — — [ — cos0 sin<^> | . (3.6b) 

sin 6 cos ( 



V2(-T + 2fi + \) 



Now we can apply the block-diagonalisation procedure (again following [5, Sec. 2.2]) to 
obtain the behaviour of eigenvalues and eigenprojections of M~ l A(e, cj))M as e — > for 
all <p. We restrict consideration to the case where A + /i ^ 0, such that Aq has two 
different eigenvalues. 

Proposition 3.1. The eigenvalues Xj(e,^) and the corresponding eigenprojections of 
A(e,(j)) have uniformly in A full asymptotic expansions as e — > 0. The main terms are 
given by 

Me,0) = T + 2X 3 + ^ + O(e% (3.7a) 

f as t + ix-\ \/2(-T + 2fX + A) L/n/ , 2 \ rt7U 
*2(e,</>) = ^ + g e + C(e), (3.7b) 

( a\ r + ^~ A \/2(-r + 2^ + A) 2 
^3(e, 0) = g e + 0(e ). (3.7c) 

Remark 3.1. The exceptional case when r = A + 2^u corresponds to isotropic media 
and is therefore excluded. In all other cases the two sheets ^(r/) = \J ^2(1]) an d ^3(77) = 
\J ^3(1]) form a double cone on the Fresnel surface S. Hence, the statement explains 
the notion of conical singularity. Note, that the linear terms are independent of 4> and 
therefore the cone is approximately a spherical cone near the conic point. 

Proof We will only shortly review the main steps. First we (l,2)-block-diagonalise 
M~ 1 A(e, 4>)M (modulo 0(e N ) for any TV" we like). The diagonaliser we are going to 
construct has the form I + eN[ ((f)) + • • • + e^ -1 ^^ ^(0) and as in [5l Sec. 2.2] its 

terms are given by recursion formulae. For we divide the off- (block-) diagonal terms 
of A\ by the difference of the corresponding diagonal entries of Aq. This gives as first 
term 

1 \ / cos (ft sin( 

M 1) W= \ -cos0 I (3.8) 



3(A + /x) 



sm ( 



and allows to cancel the off- (block-) diagonal entries of A±. We skip the further construc- 
tion and move to the next step. Since the lower 2x2 block of b-diag 1 2 A\ has distinct 
eigenvalues (namely ±1) we can now perform a diagonalisation scheme in the subspace 
corresponding to this block (modulo 0(e N )). Again we restrict ourselves to the main 
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terms. A unitary diagonaliser of the 2 x 2-block can be chosen as the unitary matrix 

,1/2(0) - I sin I cos^r j . (:t.9) 



sin 

cos 



COS 

— sin ■ 



2 2 ' 

After transforming with that matrix we apply the recursive scheme to diagonalise further. 
Note that after applying M 2 ((j)) the matrix is diagonal modulo 0(e 2 ) and therefore, 
M(I + eN^\<j)))M 2 {(j)) = MoO) + eM\{<j>) + 0(e 2 ) determines the main terms of a 
diagonaliser of the matrix A(e, (ft) and we can deduce the statements about the eigenvalue 
asymptotics. □ 

3.2.1. System formulation. Let M(e, 0) be the diagonaliser of A(e,(f>) constructed in 



Proposition 3.1 Then we consider 

^(Di + |e|P 1 / 2 (e^))M- 1 (e,0)c>(t,C) x 
V(t,e,<j>,\£\) = I (Di-|e|P 1 / 2 (e^))M- 1 (e,^)C>(t,0 



(3.10) 



with £ = |£|r/(e, (j>) and £> 1//2 (e, eft) = diag(cji(e, </>),...) the diagonal matrix containing the 
square roots Uj(e,4>) = Xj(e, <ft) of the eigenvalues of A(e, <f>). The vector V satisfies 
the first order system D t V = B(e,<f>,\£\)V with B(e,(f>,\£\) = Bi(e,<f>)\£\ + B 2 \(,\ 2 given 
by 



Bi(e, 



( wi(e, 



\-|7ai(e, 



^2(,e, ' 



ijai(e,</))\ 
\-ya 2 {e,4>) 



-uJs{e,(p 
■|7a3(e,' 



1703 (e, 




(3.11) 



/ 



and i?2 = diag(0, . . . , 0, in). The cou plin g functions Oj(e 
vector M _1 (e, cj>)r]. From Proposition 
as e — > 0. 



3.1 



are the components of the 
we know that they have asymptotic expansions 



Remark 3.2. 1. Since M _1 (e,. 
it follows that 

ai(e,4>) ■■ 
a 2 (e,0) ^ 



M|(^)(/-eAT 1 (1) ((/>))M* + C'(e 2 ) by our construction 



(3.12a) 
(3.12b) 

a 3 {e,(P) = e " 3(A + ^ ' «* — + c;«->. (3.12c) 

We know that the coupling functions vanish along three great circles through fj. We 
see that the numbering of the eigenprojections is not consistent along the circles. The 
coupling functions a 2 and 03 vanish both in the degenerate direction. 

is unitary the relation ^ ■ a 2 = 1 does not hold 



l + 0(e 2 ), 

2(h + 2X-t) . 30 , ... 2 . 
€ 3(A + ^) ' ™Y + °^ 

2(^ + 2A-r) 30 

2 v , 



2. Since we do not assume that M(e, 
for these coupling fu 
as already observed. 



for these coupling functions. However, Mq{4>) is unitary and therefore ^ ■ a 2 = 1 + 0(e) 
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3.2.2. Real and imaginary parts of eigenvalues. The coefficient matrix B(e,(j),\£\) has 
the same structure as -B(£) in Section [2j Therefore, we can conclude similar statements 
on eigenvalues and their behaviour by (a) investigating the characteristic polynomial 
and (b) block-diagonalising for small and large |£|, respectively. 

Proposition 3.2. (1) tr B(e,<p, = hc|£| 2 and det B(e, 0, |£|) = i«;|£| 2 det A(£). 

(2) B(e, (p, |£|) has purely real eigenvalues for |£| / if and only if a^ie^ (ffja^e, 4>) = 
0, i.e., e = or <f> G §Z. 

(3) 5(0, <^>, |£|) /ias i/ie rea/ eigenvalues ±^2,3(0, 0) = 3 (t + /i — A) and i/iree eigen- 
values satisfying Imu > C if |£| > c and Imy ~ |£| 2 i/ |£| < c. 

(4) T/ie quotient 

a 2 (e, 0)(zj 3 (e, 0, |g|) - Me, m 2 ) + aj(e, <t>){4,z^ 6 !gj ~ Xa(e ' f.i^ 
(4 3 (e,0,|e|)-x 2 ( e ,0)|CP)(^3( e ,0)-x3(6,0)|e| 2 ) 1 • j 

involving the hyperbolic eigenvalues ^3 of B(e, <fi, |£|) is smooth and non-vanishing 
for fixed values of |£| . 

Proof. We consider only part (2) to (4). The characteristic polynomial of -B is given by 
an expression like (2.18). If we assume that eigenvalues are purely real we can split the 
expression into real and imaginary part. We consider the imaginary part first, which 
leads to 

3 

^i 2 n^ 2 -^( e ^)i^i 2 )= o - ( 3 - 14 ) 

Therefore, real eigenvalues have to coincide with the square roots of eigenvalues of A(£). 
If we assume v 2 = Xj(e, </>)|£| 2 ) is a root of the characteristic equation, we can divide by 
the corresponding factor and obtain if e ^ (and therefore A is non-degenerate) 

a 2 (e,</>) = 0. (3.15) 

If e = the characteristic polynomial factors as 

{{v - \K\tf){v 2 - Xi|e| 2 ) - vi 2 \tf){v ~ x 2 , 3 |e| 2 ) 2 = (3.16) 

with >c\ = g(r + 2 A + Afi) and X2,3 = + fJ. — A). The first factor resembles one- 
dimensional thermo-elasticity (with r 2 = K\) and gives three roots with positive imagi- 
nary parts subject to (2.43) and (2.44). Finally (4) follows by collecting the two related 
terms in the characteristic equation of form (2.19). The imaginary part of the quotient 
is given by —k\^\ 2 in hyperbolic/degenerate directions and therefore non-zero. □ 

The quotient ( |3.13 ) may be used to determine asymptotic expansions of the hyperbolic 
eigenvalue and its imaginary part as e — > for fixed |£| and <\> g" We will follow a 
different strategy and diagonalise as e — > uniform on bounded £. 

3.2.3. Asymptotic expansion as e — > uniform in |£|. Note first, that -B(|£|, 0, (j>) is 
independent of 4> and just the system of one-dimensional thermo-elasticity (2.41) ex- 
tended by four additional diagonal entries. Since we need to understand this system 
first, we are going to recall some facts about the one-dimensional theory. As |£| becomes 
small/large we already gave asymptotic expansions of eigenvalues in Section 2.5.3 The 
bit of information which is still missing is contained in the following lemma. 
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Lemma 3.3. The coefficient matrix B(£) of the one- dimensional thermo- elastic system 
given in (2.41 ) has for £ ^ and under the natural assumptions 7, k, r > cmZy simple 



eigenvalues. 

Proof. Note that the characteristic polynomial of this matrix B(£) is given by 

^ 3 -i K |civ + r 2 iei 2 ^ + iT 2 K |ei 4 , 

which is invariant under the transform v 1— >• — v and has alternating imaginary and real 
coefficients. From that we conclude that the only possible solutions are of the form ia, 
6+ic and — o+ic for certain real a, b, c. Furthermore, from the general theory of Section[2] 
it is clear that a, c > 0. Thus, the only possibility for multiplicities to occur is if b = 0. 
Plugging in b = and multiplying the corresponding linear factors gives 



v l (\a + 2ic) 



v c 



+ 2ac) + ic 2 a. 



Comparing coefficients with the above polynomial implies that k|£| 2 
which contradicts to the positivity of all quantities involved. 



-caj (c + 2a), 
□ 



We write the coefficient matrix -B(|£|, e 
\t\- X Bi\t\^ 

where 

/ ^1 



as sum of homogeneous components in e 

(3.17) 



i?( )(|eU) + eJ B (1) (kU) + O(e 2 



£ (0) (I£U) 



17 



IjJ-2 



L02 



v-37 



-57 



17 



-W2 



-w 2 



IK 



(3.18) 



\t\J 



^ (1) (ieu) 



<*i 



-£l 



*1 



\0 



i7<52 



Sill 



3* 



1-162 30 
cos 2 



- 



Sill 



3w 



-^cos^ 










i7#2 


sin 


30 




cos 



3<t> 
2 


17^2 


sin 


30 


i7(5 2 


cos 


3<j> 
2 

/ 








(3.19) 



and lji 



3 ' 01 



2(/j+2A-r) 
3(A+/x) 



As a direct consequence of the previous lemma in combination with the asymptotic 



expansions of Section 2.5.3 we obtain 



Proposition 3.4. Assume, that A + \i 7^ and 7 2 + A + \x 7^ 0. Then the matrix 
B^{\£\,(j)) has uniformly separated eigenvalues in |£| G M, G S 1 (where ±a>2 are 0/ 
constant multiplicity two). 

Now we can apply several steps of diagonalisation based on the scheme of Sec. 2]. 
At first we apply the diagonaliser of the main part. This has only effects on the two 
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entries related to Q\ and the last row/column and deter min es the eigenvalues fo(|£|> e > <t>) 
and vf(\t;\,e,(j)) modulo e 2 . Furthermore, Proposition 3.4 allows to (1, 2, 1, 2, l)-block- 
diagonalise modulo 0(e N ), N arbitrary. 

Finally we can investigate the remaining 2 x 2-blocks and diagonalise again because 
the e- homogeneous entries ±5\e are distinct (trivially uniform in |£| and <f>). 

Proposition 3.5. Assume, that A + /i / and 7 2 + A + \i ^ 0. The eigenvalues of 
B(\£\,e,(f>) have uniformly in |£| and <f> full asymptotic expansions as e — > 0. The first 
main terms are given as 

z, (|^|,e» = 1 >o(|ei) + le|0(e 2 ), (3.20a) 
utm,e,</>) = i>fm + \i\0(e% (3.20b) 
i£^(\S\,e,<f>) = ± 1 Gh i \£\± 2 8 1 \Z\e + \Z\0(<?) (3.20c) 

where ^o(ICI) an d ^(ICI) are the eigenvalues of the one- dimensional thermo- elastic sys- 
tem with propagation speed ui\ and the signs ±i and ±2 are independent of each other. 

Remark 3.3. The statement holds true uniform in |£|. However, it is only of use as 
long as the error terms |£|e are smaller than the size of the eigenvalues. For |£| — > 
the eigenvalues of the one-dimensional thermo-elastic system behave like i^o ( I CI ) ~ I CI 2 



and v x (|C|) ~ Hence, the statement of (3.20a) is only of use if |£|e <C |£| ■> ke. if 



e <C For |£| — > 00 we know similarly ^o(|Cl) ~ |£| an d ^1 (l£l) ~ =b | ^ | , which in turn 
implies that the expansion determines the behaviour of the eigenvalues. 

This restriction is by no means a severe one; the expansion is only of interest for the 
'degenerate' eigenvalues ^ 1 3 ' ±2 (|C|, e, 0) (for which no such restriction appears). 

3.2.4. Diagonalisation for small and large To complete the picture we want to give 
some comments on expansions for small and large values of |£| under the same as- 



sumptions as in Proposition 3.5 Using the ideas from [23] we can employ the (block) 
diagonalisation scheme to separate the three non-degenerate eigenvalues from the two 
degenerate ones asymptotically and give full asymptotic expansions for them as |£| tends 
to zero or infinity. The obtained expressions coincide with the formulae from Proposi- 



tions 2.5 and 2.7 It remains to understand the behaviour of the remaining 2 x 2-blocks. 
This can be done directly by solving the characteristic polynomial as in |16l Prop. 2.7] 
or by a second diagonalisation scheme. 

We focus on the latter idea and consider the case of small |£| first. The 2 x 2-blocks 
have the form 

A(ieu*>/+( WI<U * ) _ MkU «)+°™ 2 > < 3 - 21 > 

with a function <5o(|£|> 4>) ~ e l£l- If we restrict the consideration to the zone 

Z (c) = {(|C|,e» : \£\<ce, e « 1}, (3.22) 

the remainder can be written as e|£|0(e _1 |£|) an d the standard diagonalisation scheme 
applied to the last two terms gives full asymptotic expansions in powers of e _1 |C| as 
e-^l-X), 

/o(|C|, e »±5o(|C|, e » + .... + e|C|0(e- iV |C| 7V ). (3-23) 
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A similar idea applies for large |£| in the zone 

£oo(A0 = {(l£U£> : > N, e < 1}. (3.24) 

Based on 

/ oo (|£|,e,0)I + 



<W|£|,e, 



-*»(KI,«,«) +0 « (3 ' 25) 



with a function #oo(|£|j e> 0) ~ it gives asymptotic expansions in powers of e|£| as 
e|f| -> oo. 

3.3. Cubic media, uniplanar singularities. The Fresnel surface for cubic media has 
six uniplanar singularities. Again they are equivalent and it suffices to consider the 
neig hbourhood of fj = (1, 0, 0) T G § 2 . 

We introduce polar co-ordinates near fj. Let e > and (ft G [— vr,7r). Then we set 




\fl - e 2 ( ] + ecos</> ( 1 ] + esin<£ ( ] , (3.26) 

,0; 



and use an asymptotic expansion of A{rj) as e — >• 

A(rj) = A + eA x {<t>) + e 2 A 2 ((f>) + 0(e 3 ) (3.27) 

with coefficients 

A) = diag(r, /i, /z) (3.28a) 

(cos sin <p\ 

cos (ft (3.28b) 
sin<^ y 

A 2 {(p) = {t - fi) ( cos 2 (ft J+^-^f sin 20] (3.28c) 
\ sin 2 (ft) 2 \ sin20 / 

to deduce properties of the eigenvalues and eigenprojections of A(rj) near fj. We restrict 
considerations to the case when r ^ ji. Then the following statement follows again by 
the two-step diagonalisation procedure (like in the conical case and as developed in [5], 

us). 



Proposition 3.6. Assume / 0, r ^ /j and r ^ A + 2/z. T/ien i/ie eigenvalues Xj{rf) 
and the corresponding eigenprojections have uniformly in (ft full asymptotic expansions 
as e — > 0. The main terms are given by 

x l ( v ) = T -Ce 2 + 0(e 3 ), (3.29a) 
K2 ( v ) = n + i (c + ^/c 2 - (C 2 - Z> 2 ) s in 2 (2<^ e 2 + 0(e 3 ), (3.29b) 

«S(»7) = + ^ (c " ^ 2 -(C 2 -D 2 )sin 2 (2<^ 6 2 + 0(e 3 ), (3.29c) 

where 

C= (T - /i)2 - (A + /i)2 , D = X + fj, (3.30) 

T-/i 
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Remark 3.4. This statement reflects what we mean by an uniplanar singularity. Two 
of the eigenvalues coincide up to second order. 

Proof. We follow the diagonalisation scheme. Aq is already diagonal, Ai does not contain 
(l,2)-block diagonal entries. To get expansions for the eigenvalues we have to apply two 
steps of block-diagonalisation. First we treat A\ by the aid of 



ivf)(0) 




COS ( 



sm< 



(3.31) 



such that I + eN^\(p) block-diagonalises the matrix modulo e 2 . This preserves Ao and 
= b-diagx 2 A\ and gives the new 2-homogeneous component 

r(l) _ (A + /^) 2 



A 2 + A l N[ l \ 



A l N 1 



T-H 



diag(l, — cos 



sm 



(3.32) 



The starting terms of the expansion of the first eigenvalue can be read directly from 
these matrices. For the remaining two we have to diagonalise the lower 2x2 block. This 
block has the form 

C cos 2 D sin cos < 
D sin cos C sin 2 



(3.33) 



with C, D from (3.30). Eigenvalues of this matrix are uniformly separated if the condi- 
tion 



C 2 > (C 2 -D 2 )sin 2 (2< 



i.e. C / 0, D / 



(3.34) 



is satisfied. Under this assumption the full diagonalisation scheme works through and 
the main terms can be calculated directly and give (3.29). For completeness we also give 
a unitary diagonaliser of the matrix (3.33), namely 

Moid) ' 



y / 2D 2 sin 2 (20) + 2C 2 cos 2 (20) + 2C cos 20 

/ Ccos20 + ^T -D sin 20 
D sin 20 C cos 20 + ^/1 



mi(0J m 2 ( 
-m 2 (0) mi( 



(3.35) 



where sf: = ^C 2 - (C 2 — D 2 ) sin 2 (2c; 
nuity. 



7^ f , ^ ■ Expressions are extended by conti- 

□ 



3.3.1. System form. Again we use the diagonaliser M(e, 0) of A(e,(f)) constructed in 
Proposition 3.6 to reformulate the thermo-elastic system as system of first order. For- 
mulae (3.10) and (3.11) give the corresponding representation. 
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Remark 3.5. 1. Since M _i (e 



diag(l,M 2 *(0))(/-e^ 1 (1) (0)) + O(e 2 ) in the notation 



of the proof of Proposition |3.6| it follows that the coupling functions satisfy 



ai(e,0) = l + 0(e 2 ) 

r - A - 2/i 



«2(e, < 
«3(e,' 



T-yU 
r- A-2/i 
r-yu 



(mi((j)) cos</> + rri2(4>) sin0) + C(e 2 
(mi((j)) sin0 — rri2{4>) cos0) + 0(e 2 



(3.36a) 
(3.36b) 

(3.36c) 



Since r 7^ A+2^, the function 02(0) vanishes only for (/) = k^, k £ Z, while 03(0) vanishes 
for 0= (2fc + l)J, jfe GZ. 

2. Note that in contrast to the conic situation the eigenvalues coincide to second order 
in the degenerate direction, while the coupling functions still vanish to first order (if we 
approach the degeneracy from parabolic directions). 



3.3.2. Asymptotic expansion of eigenvalues as e — > 0. We write the coefficient matrix 



i?(|£|,e, 4>) as sum of homogeneous components in e, cf. (3.17). This gives 



B<°)(|f|,. 



V-2-7 



\0 i 7 4 1} ( 



17 



17 



■37 



IK 



>J a 2 



\ 

17«2 (0) 

i 7 4 1) (0) 


i7a 2 1} (0) 

i 7 4 1) (^) 
J 



(3.37) 



(3.38) 



and £>( 2 ) (|£|,</>) has entries on the diagonal, in the last row and last column. In order 
to apply a block-diagonalisation as e — > we assume that the matrix B^°\\^\,(j)) has 
as many distinct eigenvalues as possible. This is ensured if \x 7^ r, \x 7^ r + 7 2 and 
we can (l,2,l,2,l)-block-diagonalise this matrix family. Note, that due to the structure 
of the last rows and columns, the system decouples modulo e 2 into a one-dimensional 
thermo-elastic system and one containing the elastic eigenvalues. The coupling comes 
only into play for the e 3 entries. 
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Figure 5. Cut of the Fresnel surface for hexagonal media, 772 = 0. The 
parameters are chosen as t\ = 4, T2 = 10, Ai = 2, A2 = 4, (j, = 2. The 
complete surface is obtained by rotation along the vertical axis. 

Proposition 3.7. Assume [i 7^ t , \i 7^ r+7 2 . Then the eigenvalues and eigenprojections 
of B(\£\, e, (ft) have full asymptotic expansions as e — > 0. The main terms are given by 

^o(|e|,e,0) = ^o(|C|) + |e|O(e 3 ), (3.39a) 
*f(|£|,e,0) = * 1 ± (|£|) + |£|0(e 3 ), (3.39b) 

^(ICU *) = + + 1*1^) (3.39c) 



where ^o(|£|) an d ^r~(|£|) are eigenvalues of the one- dimensional thermo- elastic system 
with parameter yfr . The signs ±1 and ±2 are independent and the parameters C and D 



are given by (3.30). 



Remark 3.6. Similar to Proposition 3.5 this statement is uniform in |£|. It will be of 
particular importance for us that the hyperbolic eigenvalues v 2 , 3 coincide up to second 
order in e with the corresponding (roots of) eigenvalues of the elastic operator. This 
will be the key observation to transfer stationary phase estimates from elastic systems 
to the thermo-elastic one. 

3.4. Hexagonal media. Finally we want to discuss the case of hexagonal media. The 



elastic operator defined by (1.4)-(1.5) is invariant under rotations in X3-direction. We 
will make use of this fact and reduce considerations to a two-dimensional situation 
corresponding to cross-sections of the Fresn el surface. 



As already pointed out in Section 2.5.4 degenerate directions are ±(0,0, 1) , which 



are uniplanar. They could be handled similar to the cubic case, but rotational invariance 
makes estimates simpler. There are further circles of degenerate directions if 

r 2 - 2n > A 2 + 2/x. (3.40) 

We exploit rotational symmetry and consider the system only in the frequency hyper- 
plane Tyi = 0. Then it is possible to express the eigenvectors rj{rf) corresponding to 
eigenvalues Xj{ri) smoothly and the thermo-elastic system can be reformulated as sys- 
tem of first order in full analogy to the general treatment in Section [2j The derived 
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asymptotic expansions for eigenvalues and the description of their behaviour transfers 
away from the degeneracy and has to be equipped with an additional description near 
these degenerate directions. 

Apart from the ^3-axis it is possible to find smooth families of eigenvectors 7^(77) of 
A{rf). This follows directly from rotational invariance combined with one-dimensional 
perturbation theory of matrices, [6]. If we assume that the frequency support of initial 
data and therefore of the solution U, 6 is conically separated from the uniplanar direc- 
tions we can follow Section [2] and rewrite as first order system in V(t, £) with coefficient 
matrix B(£) given by ( |2.9[ ) and of (l,l,5)-block structure. In what follows, we will ignore 
the scalar hyperbolic blocks and consider the remaining 5x5 matrix. 

Based on the discussions from Section |2.5.4| we know that this 5x5 block is non- 
degenerate in the sense that its 1-homogeneous part has distinct eigenvalues if (n — 
At)(T2 — /x) 7^ 0. We assume this in the sequel. But this means that the theory of 
Section [2] is applicable and gives a full discription of eigenvalues and eigenprojections 
and we are done. 

Near the uniplanar directions, i.e., on the ^3-axis, we follow the same approach as 
for cubic media. We introduce polar co-ordinates around this direction and construct 
expressions for the corresponding asymptotics. There is one major simplification, due 
to rotational invariance the construction is independent of the angular variable. 

4. Dispersive estimates 

We will show how to use the information obtained in Sections [2] and [3] to derive L p -L q 
decay estimates for solutions to thermo-elastic systems. Some of the ideas we present are 
general in the sense that they can be applied to arbitrary space dimensions, however, our 
moin focus will be the three-dimensional case and the examples considered in Section [3} 

The estimates we have in mind are micro- localised to (a) non-degenerate parabolic, 
(b) non-degenerate hyperbolic or (c) degenerate directions. The first two situations 
generalise the consideration of |16j . |22j taking also into account the estimates due to 
Sugimoto |19j , |20j , while the treatment of degenerate directions is inspired by the work 
of Liess [8], [TO]. 

4.1. Non-degenerate directions. We will consider two situations and micro-localise 
solutions to either open sets of parabolic directions or tubular neighbourhoods of compact 
parts of regular submanifolds of hyperbolic directions. 

4.1.1. Estimates in parabolic directions and for parabolic modes. Let first ip,ip G Cg (S™ _1 ) 
be supported in U with ip = 1 on supp^ and \ G C°°(M+) a cut-off function satisfying 
x(s) = for s < e and x( s ) = 1 for s > 2e. We extend both -0 and ip as 0-homogeneous 
functions to M n . Then we consider the solution to the first order system 

B t V = B(D)V, V(0,-)=i>(D)V , (4.1) 

with data V £ S(W a ;C 2n+1 ). Note, that this is well-defined and B(£) needs only to be 
defined on supp^. 
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Lemma 4.1 (Parabolic estimate). Assume that supp^ is contained in the set of para- 
bolic directions. Then the solutions to ( |4.1|) satisfy the a-priori estimates 

\\ X (\D\)^(D)V(t, -)\\ q < e- ct \\V \\ p , r , (4.2a) 

||(1 - X (\D\)MD)V(t, OH, < (1 + tyn-^\\V \\ p (4.2b) 

for all 1 < p < 2 < q < oo and mi/i Sobolev regularity r > n(l/p — 1/q). 

Sketch of proof. The proof of this estimate is straightforward from the two-dimensional 
situation considered in [16]. For small frequencies we write the solution V as sum 

V(t,x)= Y, e iMD) P»(V)Vo, (4.3) 

i/(OespecB(0 

Pt/ the corresp ondi ng eigenprojections. We know that ||-Pi,(£)ll i$ 1 an d Im^(£) ~ |£| 2 



by Proposition 2.5 Now each of the appearing terms can be estimated using the L P -LP 



boundedness of the Fourier transform (for p £ [1, 2]) and Holder inequality. Similarly, the 



representation (4.3) in combination with the bound Imspec-B(£) > C gives exponential 
decay of I? and H s norms and this combined with Sobolev embedding yields the desired 
estimate. 

For intermediate frequencies we may have to deal with multiplicities and resulting 



singularities of the spectral projections. Instead of (4.3) we use a spectral calculus 
representation which implies 

\V(t,0\ < e~ Ct ^ J IKC - 5(0)- 1 ||dC < e~ ct (4.4) 
based on the compactness of the relevant set of frequencies £ and the bound on the 



imaginary part due to Corollary 2.8 / Proposition 2.9 Here, T is a smooth curve 



encircling the family of parabolic eigenvalues for the relevant ^. □ 

If we consider hyperbolic directions we know that the parabolic eigenvalues are sepa- 
rated from the hyperbolic ones and we can use the spectral projection associated to the 
group of parabolic eigenvalues to separate them from the hyperbolic one(s). In this case 
the estimate of the above theorem is valid for the corresponding 'parabolic modes' of 
the solution. So we can restrict consideration to hyperbolic eigenvalues near hyperbolic 
directions. 

4.1.2. Treatment of non- degenerate hyperbolic directions. We consider only the for us 
interesting case when hyperbolic directions form part of a regular submanifold of S n_1 
and coupling functions vanish to first order, i.e., we assume that the corresponding 
coupling function aj : S n_1 D U — > M. satisfies 

daj(r/) 7^ when aj(r]) = 0, 77 G U. (4-5) 

This implies that Mj = {n £ IA : aj(n) = 0} is regular of dimension n — 2, the normal 
derivative d n aj(n) 7^ never vanishes and aj{rf) < e defines a tubular neighbourhood 
of Mj with a natural parameterisation. The desired dispersive estimate is related to 
geometric properties of the section °f t ne Fresnel surface lying directly over Mj, 

S [Mj) = {uj 1 ^ : n € Mj} = Sj n co Mj. (4.6) 
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Here co Mj denotes the cone over Mj. For dimensions n > 4 we have to distinguish 
between different cases, depending on whether the cross-section <S(Mj) °f the Fresnel sur- 
face satisfies a convexity assumption or not. By the latter we mean that any intersection 
of Sj with a hyperplane tangent to co Mj is convex in a neighbourhood of Srj^.y 

If this convexity assumption is satisfied (or if n = 3 and therefore dimMj = 1), we 
define the convex Sugimoto index of <S(a/,-) as maximal order of contact of <S(Af) with 
hyperplanes normal to co Mj . 

Theorem 4.2 (Hyperbolic estimate, convex case). Assume that ip is supported in a 
sufficiently small tubular neighbourhood of the regular hyperbolic submanifold Mj and 
that S^Mj) satisfies the convexity assumption. Let further jj = 7(5^)) be defined as 
above. 



Then the solutions to (4.1) satisfy the a-priori estimate 



||^(D)P^(D)^,-)|| g <(l + t)^ (|+ "" )( ^'VolU (4-7) 
for all p G (1, 2], pq = p + q and with Sobolev regularity r > n(l/p — 1/q). 

Proof. First, we outline the strategy of the proof. We split variables in the tubular neigh- 
bourhood of the regular hyperbolic submanifold Mj, one coordinate being the defining 
function aj(rj) and the other parameterising points on Mj. We have to combine a (sim- 
ple) parabolic type estimate in normal directions taking care of the imaginary part of 
the phase with stationary phase estimates for the integration along Mj. The stationary 
phase estimate is done first and follows the lines of [19], [20] along with [181 Sect. 5]. 

It is sufficient to show the estimate for t > 1. We follow the treatment of Brenner [2] 
and decompose the Fourier integral representing the corresponding hyperbolic modes of 
the solution V into dyadic pieces. For large and intermediate frequencies this amounts 
to estimate for all k G No 



l k (t) = sup 



r)=e 



e i^l(^+l?|-S(0) 



f/=-e Jf)£Mj J2 fc - 1 <|g|<2 fe + 1 

n—l—r. 



(4.t 



with the notation z = x/t, £ = |^|ry, r/ ~ (fl,fj) with f) £ Mj and fj = aj{rf). The 
amplitude Pj(£) arises from the spectral projector P v . (D) and the phase Vj(£) is complex- 
valued with Im^-(^) ~ f/ 2 uniform in £ G suppxfc and k 6 No- 

If z + V£i^-(£) 0, £/|£| € Mj or if z is not near a direction from Mj, the principle 
of non-stationary phase implies and gives a rapid decay. It suffices to restrict to z 
corresponding to stationary points. We use the method of stationary phase to estimate 
the integral over Mj, this can be done uniformly over £ and f], provided e is chosen small 
enough and yields an estimate of the form 



/ . . . dr) 



< Ct ii r ii e"^ 2 * (4.9) 



uniform in k and \fj\ < e. In order to obtain this estimate we apply Ruzhansky's multi- 
dimensional van der Corput lemma, [T7] , based on the uniformity of the Sugimoto index 
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j(Sj n co {r/ : aj(r]) = fj, n ~ 77}) for small r] and the uniform bounds on the appearing 
amplitude. Similar to [16] the imaginary part of the phase can be incorporated in the 
estimate for the amplitude. Integration over fj yields a further decay of t _1//2 , while 
integrating over £ and using |£| ~ 2 fc yields 

l k (t) <Ct 2 ^ 2 1 Ti j . (4.10) 
Hence, we need r > n — (compared to the elasticity or wave equation with r > 

n — to apply Brenner's argument and obtain the desired estimate for the high 

frequency part. The required regularity follows from using Sobolev embedding for small 
t. 

The treatment of small frequencies is somewhat simpler. We do not apply a dyadic 
decomposition, but still have to use a stationary phase argument along Mj combined 
with the behaviour of the imaginary part of the phase away from it, 



T(t) = sup 







fjeMj J\£\<i 

n-2 f fV= e 9 „ o n-2 _1 

< Ct H / / e~ cfl m ^ d|f | < Ct ^ 

J\e\<lJfj=-e 



□ 



Without proof we comment on the non-convex situation. If the convexity assumption 
is violated we have to replace the convex Sugimoto index by a corresponding non-convex 
one 7o(«S(m,-))- This is defined as the maximum over the minimal contact orders of S(Mj) 
with hyperplanes normal to the cone co Mj, the maximum taken over all points of <S(jw-). 
The price we have to pay for non-convexity is a loss of decay. 

Theorem 4.3 (Hyperbolic estimate, non-convex case). Assume that ip is supported in 
a sufficiently small tubular neighbourhood of the regular hyperbolic submanifold Mj and 
that S^Mj) does not satisfy the convexity assumption. Let further jj = jo{S(Mj)) be the 
non-convex Sugimoto index. 



Then the solutions to (4.1) satisfy the a-priori estimate 



D)P Vj (B)V(t,-)\\ q <(l + t) (|+ ^ )( - «VolLr (4.11) 

for all p G (1, 2], pq = p + q and with Sobolev regularity r > n(l/p — 1/q). 

4.1.3. Application to cubic and hexagonal media. Because of its importance later on we 
remark that in our applications to three-dimensional thermo-elasticity the manifolds Mj 
are parts of circles on § 2 , i.e. can be seen as intersections of §> 2 with a cone. So we 
have to look at the corresponding sections of the Fresnel surface. In this case 7^ is just 
the maximal order of tangency between the curve 5(j^.)and its tangent lines. If the 
curvature of this curve is nowhere vanishing, then 7,- = 2. Furthermore, algebraicity of 
S of order 6 implies that the highest order of contact is 6 and therefore jj G {2, . . . , 6} 
is the admissible range of these indices. 

For cubic media there are two types of regular hyperbolic submanifolds. One is up 
to symmetry given by the circle 773 = on S 2 and the corresponding eigenvalue is 



28 



JENS WIRTH 



equal to [i. Thus the section of the Fresnel surface is just a circle and therefore its 
curvature is nowhere vanishing. Similarly, for intersections of the Fresnel surface with 
the plane m = ^3 we obtain the hyperbolic eigenvalue >c = "n\{j — X) + rjffjL. It is a simple 
calculation_]to show that the curvature of the corresponding section of the Fresnel surface 
is nowhere vanishing as soon as A / r and \i 7^ 0. Hence, jj = 2 in both cases. 

For hexagonal media regular hyperbolic submanifolds correspond to circles on the 
Fresnel surface. Again, 7,- = 2. 

4.2. Cubic media in 3D. We want to discuss the derivation for estimates near degen- 
erate directions by the example of cubic media in three-dimensional space and combine 



them with the general estimates from Section 4.1 



4.2.1. Conic points. The following statement resembles [H Thm. 1.5]. In [10, Sect. 3] a 
stronger decay rate is obtained for some conic degenerations, but they require a suffi- 
ciently bent cone and we can not guarantee that in our case. 

Theorem 4.4 (Conic degeneration). Assume U\, U2 and 9q are micro-locally supported 
in a sufficiently small conical neighbourhood of a conically degenerate point on S 2 . Then 
the corresponding solution to the thermo- elastic system for cubic media satisfies the a- 
priori estimate 

\\y/A{p)U(t, •), U t (t, -),e(t, -)\\ q < (1 + ty^-^\\ y r A(D)U 1 , U 2 ,8 \\ P ,r (4.12) 
for p € (1, 2], pq = p + q and r > 3(1/ p — 1/q). 

Proof The main idea is that the proof of [7{ uses polar co-ordinates around the singu- 
larities of the Fresnel surface similar to our treatment in Section [3| Stationary phase 
arguments are applied in tangential direction and are uniform for small radii, while the 
final estimate follows after integration over the remaining variables. 

It suffices to prove the statement for t > 1, the small time estimate is a direct conse- 
quence of Sobolev embedding theorem in combination with the obvious energy estimate. 
Similar to the hyperbolic estimate discussed before, we apply a dyadic decomposition of 
frequency space (localised to a small conic neighbourhood of the degenerate direction). 
The estimate for single dyadic components follows |7] resp. [8, Thm. 1.5]; the only thing 
we have to check is that the necessary assumptions are satisfied uniform with respect to 
|£| and k £ N. We consider 



/■e /•2-7T 

Z k (t) = sup 



e 

JO 72 fc - 1 <|^|<2 fc + 1 



iti€K*-»?+[erS(i£i,e,*)) 



x ^(k|,e»Xfc(iei)ier" r d|e|d06de 



, (4.13) 



where rj £ E 2 denotes the point with polar co-ordinates (e, 4>) near the conic degenerate 
direction and £ = \^\m The amplitude pj(|£|,e, (j>) arises from the spectral projector 
(given in terms of the diagonaliser) constructed in the blown-up polar co-ordinates and 



^Parametrising by the angle, the hyperbolic eigenvalue is given by *e(<f>) = fi + T A 2 2m 
remains to check that 9? \J x(<j>) + \J x(</>) 7^ 0, see [22] for such a calculation. 
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Xk(0 corresponds to the dyadic decomposition. The complex phase z^-(|£|, e, <fi) is de- 
scribed in Proposition 3.5 Its imaginary part is non- negative and vanishes to second 
order in e = 
erate point 



as well as for three hyperbolic manifolds emanating from the conic degen- 
Again we may treat this imaginary part as part of the amplitude and apply 
stationary phase estimates for the integral with respect to eft. As the approximation of 
the phase modulo 0(e 2 ) is independent of (ft and uniform in |£| this yields 

'<H|£|§-*- e t (4.14) 



uniform in |£|, k and < e < e. There is no further benefit from the imaginary part (as 
there can not be a lower bound with respect to e) and integrating with respect to |£| 
and e concludes the estimate for Xj^(t). Similarly, we estimate the small frequency part 

f2n 



X(t) = sup 



J 



KI<1 



3 i ^l(^+^- 1 ^(l«l^))p i (|e|,e^) X (|^|)|C| 2 d|e|d0ede 



< cr 



such that Brenner's method again yields the desired decay estimate. 



(4.15) 

□ 



4.2.2. Uniplanar points. The treatment of uniplanar degeneracies follows [8]. We have 
to make one further additional assumption related to the shape of certain curves on the 
Fresnel surface near the degenerate point. To be precise, we either require that 

Q fl S n IT has non- vanishing curvature (4.16) 

for fi C M n an open neighbourhood of the uniplanarly degenerate point and for any 
plane n sufficiently close and parallel to the common tangent plane at the unode. This 



condition is equivalent to the technical assumption (1.12) made in |10| If (4.16) is vio- 
lated, we need to consider Sugimoto indices ^ u = 7(f2 D S fl II C II), i.e., contact orders 
of these planar curves with its tangent planes combined with a uniformity assumption. 
Under assumption (4.16) the index is given by j u = 2. 



For cubic media we have to use the statement of Proposition 3.6 to determine the 



index j u . Using the notation of (3.30), it suffices to calculate the indices of the indicator 
curves determined by e 2 (fi + C ± \/C 2 cos 2 (20) + D 2 sin 2 (20)) = 1. This yields 

7«G {2,3,4} (4.17) 

In the nearly isotropic case we have 7 n = 2, away from it 7^ = 3. Both are generic, while 
the borderline case with j u = 4 is not. The asymptotic construction of the eigenvalues 



and eigenprojections near the uniplanarly degenerate point of Proposition |3 . 7| yields that 
the assumption is satisfied uniformly for the phase functions appearing in all dyadic 
components of the operator. 

Theorem 4.5 (Uniplanar degeneration). Assume U\, U2 and Oq are micro-locally sup- 
ported in a sufficiently small conical neighbourhood of a uniplanarly degenerate point on 
§> 2 . Let further be the index of the uniplanar point. Then the corresponding solution 
to the thermo- elastic system for cubic media satisfies the a-priori estimate 

\\y/A(p)U(t, •), U t (t, -),9(t, OH, < (1 + t)-^ + ^p-^\\y/A^)U lt U 2 M\p,r 
for p G (1, 2], pq = p + q and r > 3(1 /p — 1/q). 



(4.18) 
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Sketch of proof. We will sketch the major differences to the treatment of conic degenera- 
cies. We will again use polar co-ordinates and estimate corresponding dyadic components 



(4.13), where now e, 4>) is determined by Proposition 3.7. The imaginary part of 

Vj(\^\,e,(j)) vanishes to third order and is of no benefit, while the real part coincides to 
third order with the corresponding elastic eigenvalue. This allows to use estimates from 
[8] and [101 Sect. 4], the main difference to the previous situation is that we now use 
stationary phase estimates for both, the angular and the radial integral. The proof itself 
then coincides with the corresponding proof for cubic elasticity, cf. [H] . 

Using a change of variables the integral is written in the new variables cjj(r/)|^| (i.e., 
roughly Heuj) and rj/ujj(r]) £ Sj. In this form the phase splits and the crucial estimate 
is just a Fourier transform of a density carried by the sheet of the Fresnel surface (with 
possible singularity in the unode) . This is calculated by introducing distorted polar co- 
ordinates on the surface. As level sets we use cuts of the surface by planes parallel to the 
common tangent plane. Then we will at first apply the method of stationary phase to the 
radial variable in these co-ordinates. These stationary points are non-degenerate and we 
use the obtained first terms in the asymptotics for a second stationary phase argument 



in the angular variables. The condition (4.16) would imply again that stationary points 



are non-degenerate and we are done, while if (4.16) is violated we use the lemma of van 
der Corput instead to prove the estimate. □ 

4.2.3. Collecting the estimates. It remains to collect all the estimates into a final state- 
ment for cubic media. Parabolic directions are treated by Lemma 4.1 hyperbolic man- 
ifolds away from degenerate points are covered by Theorem 4.2 The remaining 24 
degenerate directions fall into either of the previously discussed categories and estimates 
follow from Theorem 4.4 and 4.5 The resulting estimates are collected in Table [TJ 





small frequencies 


large frequencies 


parabolic directions 


(l+t)-2 


e -ct 


hyperbolic directions 




(1 + t)" 1 


conic degeneracies 




(l + i)-3 


uniplanar degeneracies 


7 G {2,3,4} 


(i + t)~M 

7 G {2,3,4} 



Table 1. Contributions to the dispersive decay rate for cubic media. 
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Corollary 4.6 (Cubic decay rates). Cubic media in three space dimensions satisfy the 
dispersive type estimate 

\\y/Afp)U(t, •), U t (t, -),8(t, Olk^K") < (1 + tr* ( *-«>|| s/MpjUx, U 2 , 9 \\ p ,r (4.19) 

for all data U x G WP> r+1 (R 3 ; C 3 ), f/ 2 G W^' r (IR 3 ;C 3 ) and 6> G W^ r (IR 3 ) ; provided p G 
(1, 2], = p + g and r > 3(l/p — 1/g). 

Decay rates improve if the Fourier transform of the initial data vanishes in the conically 
degenerate directions. This could be achieved by posing particular symmetry conditions. 



4.3. Hexagonal media. The treatment of hexagonal media is somewhat simpler. The 
uniplanar degenerations trivially satisfy the assumption (4.16) and therefore yield the 
decay rates specified by the above theorem. The additionally appearing manifolds of de- 
generate directions are trivially resolved as there are smooth families of eigenprojections 
associated to both eigenvalues (as we stay away from the uniplanar points) and we can 
therefore treat the modes separately. 

One of them is hyperbolic for all directions, we refer to it as the genuine hyperbolic 
mode. The sheet of the Fresnel surface corresponding to this mode, i.e., to the eigenvalue 
x ( r ?) = n ~2 Kl (vf + V2) + wl i s eas ily seen to be strictly convex for all choices of the 
parameter and gives therefore t . The proof is similar to that for the wave equation, 
see 0. 

The parabolic modes away from the degenerate hyperbolic directions are treated as 
before, while the remaining degenerate hyperbolic manifold is treated by the estimate 
of Theorem |4.2| with 7 = 2 due to rotational invariance. The resulting estimates are 
collected in Table [2 





small frequencies 


large frequencies 


genuine hyperbolic mode 


(1 + t)- 1 


(1 + t)-' 


parabolic modes 


(1 + t)- 3 / 2 


e -ct 


hyperbolic directions 


(1 + t)" 1 


a+r 1 


uniplanar degeneracies 


(1 + t)" 1 


a+r 1 



Table 2. Contributions to the dispersive decay rate for hexagonal media. 
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Corollary 4.7 (Hexagonal decay rates). Cubic media in three space dimensions satisfy 
the dispersive type estimate 

\\y/A(P)U(t, •), U t (t, -),9(t, -)||l«(r») < (1 + ty^W^MB)^, U 2 , 6o\\ p , r (4.20) 

for all data Ui £ W p ' r+1 (R 3 ; C 3 ), U 2 € W p > r (R 3 ;C 3 ) and 6 G W™ r (R 3 ), provided p G 
(1, 2], = p + g and r > 3(1/ p — 1/q). 
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